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Noncompact S0(1, A^) sigma-models are studied in terms of their large ex- 
pansion in a lattice formulation in dimensions d > 2. Explicit results for the 
spin and current two-point functions as well as for the Binder cumulant are 
presented to next to leading order on a finite lattice. The dynamically gen- 
erated gap is negative and serves as a coupling-dependent infrared regulator 
which vanishes in the limit of infinite lattice size. The cancellation of infrared 
divergences in invariant correlation functions in this limit is nontrivial and 
is in (i = 2 demonstrated by explicit computation for the above quantities. 
For the Binder cumulant the thermodynamic limit is finite and is given by 
2/{N+l) in the order considered. Monte Carlo simulations suggest that the 
remainder is small or zero. The potential implications for "criticality" and 
"triviality" of the theories in the S0(1, A^) invariant sector are discussed. 
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1. Introduction 



In quantum field theories witli nonabelian symmetries and dynamical mass generation 
the large N expansion often provides a qualitatively correct and quantitatively rea- 
sonable description of the physics of the systems. Specifically in sigma-models with a 
compact global symmetry group the expansion is known to be an asymptotic expansion 
PP and when slightly ad hoc applied to low orders at fixed small N sometimes gives 
surprisingly accurate results, see e.g. [21 for the renormalized coupling. In a lattice 
formulation one starts off on a finite lattice, the associated 'finite volume' mass gap then 
is uniformly bounded away from zero, and in the large N series for invariant correlation 
functions the limit of infinite lattice size (also called the thermodynamic limit) can safely 
be taken termwise. 

A study of the large N expansion in sigma-models having a noncompact S0(1, A^) in- 
ternal symmetry group has been initiated in [J| fS]. A large value of A^ in this case is 
also physically relevant for the granular limit of random hamiltonians describing disor- 
dered electrons with A^ orbitals per site [6]. In a lattice formulation of the S0(1,A^) 
sigma-models again a gap is dynamically generated in the large A^ expansion, which 
is however negative and vanishes as the size of the lattice goes to infinity. Effectively 
the gap now acts as a subtle, coupling- dependent, infrared regulator and the techni- 
cal problem consists in studying the 'coordinated' limit V ^ oo of lattice sums of the 
form k„ fv{ki, . . . , kn), where fv carries an explicit ^-dependence via the gap. 

The sums associated with individual Feynman diagrams of the large A^ expansion will 
typically diverge in the limit. The issue whether or not in the combinations entering 
invariant correlation functions the infrared divergences cancel is analogous to the one 
encountered in the perturbation theory of compact sigma-models |71,[8j and is the subject 
of the present paper. Since this issue is most critical in the two-dimensional systems we 
examine the limit specifically in this case, although our finite volume results are valid 
in all dimensions d > 2. We compute a number of physically interesting quantities to 
leading and subleading order and show that they indeed do have a well-defined ther- 
modynamic limit. Concretely we consider the spin two-point function, the two-point 
function of the Noether current, and the Binder cumulant. 

The Binder cumulant U is defined in terms of the zero momentum limit of the connected 
four-point function. In massive scalar field theories it serves to define an intrinsic measure 
of the interaction strength and has been used to explore "triviality" issues. In a massless 
theory, like the systems considered here, there is no obvious reason why U should have 
a finite thermodynamic limit. Somewhat surprisingly we find that U does have a finite 
and nonzero limit to leading and subleading order, which is moreover independent of A 
and given by 2/(A^+l). Supported also by Monte-Carlo simulations we conjecture that 
the infinite volume limit of the exact U is also very close to 2/{N + l). The potential 
implications for "criticality" and "triviality" in the S0(1,A^) invariant sector of the 
theory will be discussed in the conclusions. 

The rest of the article is organized as follows. In the next section we review a result from 
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a previous paper [5] which on a finite lattice allows one to do large computations in 
the simpler compact models and then transfer the results to the noncompact ones via 
a "large N correspondence" . This correspondence has an interesting interplay with the 
Schwinger-Dyson equations and instead of going through the (fairly routine) diagram- 
matic computations we merely present the results as solutions of the large expanded 
Schwinger-Dyson equations with the correct 'initial' data. Expressions for the two- and 
four point functions to leading and subleading order are given (valid on a finite lattice 
in all dimensions d >2), from which also the two-point function of the Noether current 
and the Binder cumulant can be obtained in the same order. The thermodynamic limit 
in (i = 2 of the local quantities and of the Binder cumulant are studied in Sections 4 and 
5, respectively. 



2. Large N expansions of compact and noncompact models 

In compact sigma-models the large N expansion is a saddle point expansion based on a 
generating functional obtained by 'dualizing' the spins, i.e. by imposing the constraint 
via a Lagrange multiplier field and performing the Gaussians. The counterpart of this 
duality transformation is somewhat ill-defined in the noncompact models. The large 
expansion can nevertheless be justified and on a finite lattice the expansion coefficients 
for invariant correlation functions can be inferred from those in the compact model [3]. 
This "large correspondence" allows one to do computations in the compact model, 
where no gauge gauge-fixing is required, and the familiar framework can be used. Here 
we briefiy summarize the correspondence and present explicit results for two and four- 
point functions in Section 3. The results of Sections 2 and 3 are valid in all dimensions 
d>2. 



2.1 Definitions 

Here we recall the notation and the definitions for the invariant correlation functions 
considered and their generating functionals. We consider the SO(A^ + 1) spherical and 
the S0(1,A^) hyperbolic sigma-models in two dimensions with standard lattice action, 
defined on a hypercubic lattice A C Z*^ of volume = |A| = L'^. The dynamical variables 
("spins") will be denoted by n^, x G A, a = 0,...,N, in both cases, and periodic 
boundary conditions are assumed throughout rix+Lfi = n^. The constraint is n ■ n = 1 in 
both cases, but with different 'dot' products; namely a-h := oPlP +a^h^ + . . .+a^b^ in the 
compact model, and a - b := aPlP — a^b^ — ... — a^6^ in the noncompact model. Clearly 
= {ne M^+^ I n ■ = 1} is the A^-sphere and H'^ = {n e M^-^ \ n ■ n = 1, > 0} is 
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the upper half of the two-sheeted A^-dimensional hyperboloid. The lattice actions are 



S± = • ^^+A - 1) = ^ ri^ • (An)^ > , (2.1) 



where the upper sign refers to the compact model and the lower sign to the noncompact 
model. The laplacian is A^y — — '^^j^i'^S^^y — 5x,y+ix ~ ^x,y-jx\-i ^ usual. We write 

dO+(n) = d^+^n(5(n-n-l), 

dQ_(n) = 2d^+^n5(n -71-1)^(71°), (2.2) 

for the invariant measure on and , respectively. Further 5±{n, n') is the invariant 
point measure on 5"^, , and v) — (1,0, ...,0). Note that the measure df^+(n) is 
normalized while has infinite volume. 

In the compact model we consider the generating functional, 

exp W+[H]=MjY[ dfi+(n,) cxp { - 5+ + ^ ^ H,y{n, ■ - 1) } , (2.3) 



x,y 



where Hxy > is a source field and the normalization M is such that H^[0] = 0. For the 
noncompact model we consider the generating functional 

e^^W_[H] = AA/"j]dl]_(n,)(5_(n,3,nT) exp{-5_ + ^^i/,^(n,-n,-l)},(2.4) 

X x,y 

where now H^y < sources give damping exponentials, and one spin at site Xq is fixed 
in order to make the generating functional well defined. 

Connected 2r point functions are defined by 

W±[H] = W±^rixi, yi; . . . ; Xr, Vr) H^,y, ■ ■ ■ H^^Vr > 

W±^r{.Xi,yi;...;Xr,yr)-=K^y^---h^ryrW±[H] , Ky ■= YT]— ' (2-5) 

H=0 OHxy 

In particular W±^i{x,y) = {n^c ■ ny)± - 1, PF±,2(a;i, 2:2, 2/2) := {uxi ■ riy^n^^ ■ ny^)± - 
{uxi • ny^)±{nx2 ■ 'nyr^)±^ where ( )± are the functional averages with respect to J\f~^e^^^. 
Note that VF±,r(- . . ; x, x; . . .) = 0. 
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2.2 The l/N expansion 



The goal in the following is to construct these invariant correlation functions in a large 
asymptotic expansion. That is, A := {N+1)/ (3 is kept fixed and the coefficient functions 
W±^r in 

A^ °° 1 

W^±,r(a:i, . . . , y.) = ^^^^y., E (iv + i)s ^£^^^^' • • • ' ' ^2.6) 

are sought, with the understanding that the right hand side of (12.61) provides a valid 
asymptotic expansion of the exact W^-, initially on a finite lattice. 

is) 

The diagrammatic algorithm for the computation of the coefficient functions W\_ ^ is 
rather straightforward in the compact model, see e.g. [9|. From |lj it is also known to 
provide a valid asymptotic expansion. Direct computation of the functions w'f \ in the 
noncompact model is also possible ^ , although due to the gauge fixing the computations 
are considerably more tedious than in the compact model. In [10] it will be shown that 
this algorithm also provides a valid asymptotic expansion (12. 6 p for the W-^r- 

One of the advantages of a lattice formulation for these systems is that there is an 
exact correspondence [5] between the functions w'^\ in the noncompact model and their 

counterparts W\ ^ in the compact model, valid on a finite lattice in all dimensions d > 2: 

(a) The coefficient functions PV±,r are translation invariant and can be expressed in 
terms of D±{x) = D{x)\uj->uj±, with D{x) the free propagator of squared mass u and 
with uj±{X,V) the solutions of the gap equations XD{0) = ±1 discussed further in 
Subsection 3.1. 

(b) For all r > 1, s > 0, there exists unique functionals X^'\d]{X) oi D such that W^'l = 

xi^\D+]{X) are the coefficients in the compact model and wj:^l = (— l)''Xr'*'*[D_](— A) 
are the coefficients in the noncompact model. 

As a consequence the computations only have to be done in the compact model, and 
the result in the noncompact model can be obtained via (b). In the next section we will 
compute a subset of correlation functions, but instead of doing the computation using 
the 1/A^ Feynman rules we present the results and verify that they solve the associated 
Schwinger-Dyson equations. 
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3. Schwinger-Dyson Equations 



In the compact model the Schwinger-Dyson equations for the functions ( 12. Sp have, to our 
knowledge, first been formulated by M. Liischer [llj. The derivation is readily extended 
to noncompact models and is reproduced in [5j. Here we just record the basic equations: 



+ ^ H^, h^yW± - h,^W± - KyW± - KzKyW± - {h^^W±){h^yW±) 

-N{l-6,y)iKyW± + l)=0, 



(3.1) 



where the upper sign corresponds to the compact model and the lower sign to the 
noncompact model. In terms of the multi-point functions (12. 5p these equations amount to 
an infinite coupled system of nonlinear partial differential equations. As such boundary 
conditions have to be specified; without them even the exact equations (13.11) do not 
determine their solution uniquely, see ^ for a counter example. However, since (13. ip 
does not contain a closed equation for any of the Wr, it is difficult to impose such 
boundary conditions in practice. 

In contrast, the large ansatz (12. 6p effectively converts the Wr equations into a hierarchy 
which can be solved recursively and where 'initial' conditions can be specified. The 
recursion pattern for the r + s > 1, functions is given in Fig. 1. To compute a 

given coefficient all quantities having arrows pointing towards it are needed. 



i i 
i 



Fig. 1: Recursion pattern for the solution of the large expanded SD equations. 

The first few equations for the W^^ are spelled out in Appendix A; a closed formula for 
the generic equation can also be given and used to show the recursion pattern in Fig. 1 
by induction. The key assumption in our use of the large A^ expanded Schwinger-Dyson 
equations will be that at each recursion step in Fig. 1 there exists a solution and that 
the solution is unique. This assumption implies that, once W^'^^ has been specified, there 
will be an infinite sequence of functions r + s > 1, uniquely associated with it, 

which in turn determine the series (12. 6p uniquely for each W,.. The choice of W^^^ is 
ultimately determined by the physics problem one seeks to study; different choices are 
possible [5] for the same (initial) equation flA.ip . 
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The existence and uniqueness of a solution at each recursion step of Fig. 1 is presumably 
difficult to establish directly from the equations. In terms of the underlying discretized 
functional integral (which solves the exact equation (13.11) by construction) existence and 
uniqueness of the solution to the recursive equations amount to the existence of a well- 
defined asymptotic expansion of the form (12.61) . For the generating functionals 
and the latter is guaranteed by the results of [1] and [10], respectively. This 

provides an indirect justification of the assumption stated in the preceding paragraph. 

It also provides the rationale for the procedure adopted in Subsections 3.2 and 3.3: we 
present the expressions for and to leading and subleading order (s = 0, 1) and 
claim that they solve the Schwinger-Dyson equations with the correct initial conditions, 
Wi^l, respectively. The required equations are tabulated in Appendix A. The verification 

that the W^l presented indeed solve these equations is straightforward and is omitted. 
It is however far shorter than the diagrammatic computation in either the compact or 
the noncompact model. 

Though in this paper we shall be concerned with the large N expansion exclusively, 
let us add that the Schwinger-Dyson equations (13. ip can also be subjected to a pertur- 
bative expansion, i.e. the ansatz (12. 6p is replaced with one in terms of powers of 1//3. 
The recursive pattern determining the coefficients is similar to that in Fig. 1, but the 
differential part of the equations now involves a linear differential operator with constant 
coefficients. This has been used by Liischer [11] to show directly from the equations that 
each recursion step has at most one solution. The existence of a solution of course follows 
from the diagrammatic algorithm as described (on the lattice) by Hasenfratz Once 
the results are known the perturbatively expanded Schwinger-Dyson equations provide 
an efficient way to verify them. 



3.1 2— and 4— functions to leading order 

The leading order two-point functions provide the starting point for the recursion. 
They have to solve Eq. flA.ll) but in order to pick a specific solution further information 
has to be added. The appropriate solutions turn out to be given by 

Wi%, y) = ±D±{x - y) - , (3.2) 

with 

^±(^) = T7E^^^' (3.3) 
V Ep + iu± 

where the sum goes over p = ^{rii, . . . , rid) , = 0, 1, .., L — 1 and Ep = X])!=i Pfi with 
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= 2 sin Further a;± are the particular solutions to the 'gap equations' 

±A-i = D±(0) = i5^-^, (3.4) 

V Ep + uj± 



obeying a;+ > in the compact model and > tu- > —3;^ sin^ vr/L in the noncompact 
model |5j. The rationale for the choice of these solutions of flA.ip is that their properties 
are necessary for the stability of the expansions, see pLj for the compact and [lOj for the 
noncompact model. All W^^^j., r + s > 1, are then in principle uniquely determined by 

the wf\^ and we shall simply present the solutions of the associated Schwinger-Dyson 
equations. 

The solution of (1A.3P for the 4-point function in leading order is 

W^l{xi,yi; X2, y2) = D±{xi - X2)D±{yi - ^2) + ^±(2:1 - y2)D±{yi - X2) 

-2 ^ - u)D^{y^ - m) A±(m - v)D^{y - X2)D^{v - 1/2) , (3.5) 

u,v 

where A±{u — v) is defined by 

J2D±{x-ufA±{u-v) =6^^, (3.6) 

u 

so that 

1 



(A:) 



Obviously in the compact model n+(fc) > for all k. In contrast in the noncompact 
model one has n_(/c) < ,k ^ and n_(0) > 0, a fact of importance for the validity of 
the large N expansion [S]. A±{u) is the propagator of the auxiliary field in the functional 
treatment of the 1/N expansion, and correspondingly the last term in fl3.5p corresponds 
to a tree diagram with an intermediate auxiliary field propagator. 



3.2 2— and 4— point functions to next-to-leading order 

The solution of flA.2l) for the 2-point function in next-to-leading order is 
Wi%,y) = -2q±-^wi%,y) 

T2 J2 D±{x - u)D±{u - v)Ai{u - v)D±{v - y) , (3. 
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where 

<1± = 1^ J]/^±(w)A±(w-t;)Z}±(w-t;)Z}±(t;) 

n±(o)\/2 E E + ^^)2 (^^^^ ^ n±(g) ■ ^^-^^ 



and the partial derivative means ±A^n±(0)^ (at fixed volume). In particular 

:= yD±{w-x)D±{w-y) = ^Y— (3.10) 

Note n±(0) = ^±,2(0). The first term on the rhs of (13. 8p corresponds to the tadpole 
diagram and the second term to the non-trivial self-energy diagram. 

Finally the solution for the next-to-leading 4-point function is given by 

d 

W^l{xi,yi;x2,y2) = -'^q±-Q^W^l{xi,yi;x2,y2) 

-2 J2 W^lixi,yi;u,v)D±{u-w)A±{u-w)Dg\v-z)W^l{x2,y2;w,z) 

u,v,w,z 

+2 J2 W^l{xi,yi,u,v)D±{u-v)A±{u-w)A±{v-z)D±{w-z)W^l{x2,y2-,w,z) 

u,v,w,z 

- E z/i; u, ^) A±(m - v)wi%2, y2; u, v) . (3.11) 

u,v 

Again the separate contributions in (13.111) become more apparent when drawn as corre- 
sponding Feynman diagrams. 

3.3 Two— point current correlation function 

In both models the Noether currents are given by 

Jf{x) = (3 [nld.n', - nld,nl\ , < a,b < N . (3.12) 
The invariant two-point function of these currents is 



+ ! 



J.M^, y) - { . (3-13) 
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It obeys the Ward identity 

d;j±M^, y) = ±2NPE± {6^,,y - 6,,y+,) , (3.14) 

where 9*/(x) = f{x) — f{x — fi) and in the noncompact case (13.141) holds for x xq 
only. Further 

E± = {n^.-n^+o)±, (3.15) 

is independent of x because of translation invariance. One way of obtaining (13.141) is 
by specializing the ' pre- Schwinger- Dyson' equation Eq. (3.39) of [5J to C = J'lf'iy) and 
using the completeness relations Eq. (3.38) in [S]. 

The two-point function can be expressed in terms of the 2- and 4-point functions of the 
spins according to 

J±,t,u{x, y) = 2/5^ W±^2{x, y;x + fi,y + u) ~ W±^2ix, y + u;x + fi,y) 

+W±,i{x, y)W±,i{x + fx,y + u)- W±,i{x, y + u)W±,i{x + fi, y) 

+W±,i{x, y) + W±,i{x + fi,y + v)~ W^±,i(x, y + v) - W±,i{x + /i, y)] . (3.16) 

The current correlation function has accordingly a 1/A^ expansion of the form: 

s>0 ^ ' 



In the lowest order we have 



J±A^. y) = pS(a;, y) + W^i^Kx + /i, y + i>) - W'^^^ix, y + u) - W^^l{x + /i, y) 

+1^1°; (x, i/)iy£}(x + A, y + A) - l^£;(x, y + z>)iy£l(x + /i, y) . (3.18) 
Inserting the solution (13.21) one gets 

ji°),^(x,y) = D±{x - y)D±{x + fi- y-O) - D±{x -y- v)D±{x + fi- y) , (3.19) 
the Fourier transform of which is 

X 

fi \ 2 sin(p + g/2)^sin(p + g/2)^ 
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It is seen to satisfy 



(1 -e-^"?-^) D±(z>) 



(3.21) 



which is the Ward identity (I3.14p to lowest order 1/A^. 
In next order we have 



J±iu{X' y) = J±Ui^^ y) + ^±^2(3;, y;x + fi,y + 0) - W^^'^i^, y + i';x + fi,y) 
+iy£i(x, y) \w^\{x + fi,y + i)) + \-^] + W^\{x + fi,y + i>) \w^l{x, y) + 

-W^S(x,y + z>) [w^i°i(x + /i,i/) + A^^] -W^^,{x + fi,y) y + z>) + A'^ 

Using the solutions (13. 2p . (13. 5p . (13. Sp . this becomes 



r(o) 



.(0) 



(0), 



(3.22) 



J±!t,v{x^ y) = 2g± ^±,2(3; - y)D±{x + fl-y -O) + D±{x - y)D±^2{x + jl - y - u) 
-D±^2{x + A - y)D±{x -y-v)- D±{x + fi- y)D±^2{x -y-v) 
-2 ^ /\±{w)D±{w) £"±,2(1/ -x + w)D±{x -y + jl-u) 

w 

+D±,2{y - X + v - fi + w)D±{x -y) - D±^2{y - x - fi + w)D±{x -y-v) 
-D±^2{y - X + + w)D±{x - y + p,) 

-2 ^ £)±(x - u)D±{y - u)A±{u - v)D±{x + fi - v)D±{y + - v) 

u,v 

+2 ^ D±(x - u)D±{y + 0- u)A±{u - v)D±{x + fi - v)D±{y - v) . (3.23) 

u,v 

Its Fourier transform is simpler in form: 

J±lu{q) = - exp ('-[q^ - q,]] [X±,i.^,(g) + X±,2;M.(g) + X^^^-Mq)] , (3.24) 
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with 

8q± ^ sin(p + sm{p + 



V ^ (£;p + ^±)2(£'p+, + ^±) 

sm(pi + 5/2)^sm(pi + 

Pl,P2 



X (a) = — V sm{pi + q/ 2)^sm[pi + q/ 2)^ 



It can be checked to satisfy the Ward identity in the next order in the large expansion. 



3.4 The Binder cumulant 

In scalar field theories with a mass gap a renormalized 4-point coupling is defined in 
terms of the Binder cumulant 

^ - ^ + iVTl " " ~ 5] {nx^-ny,n^,-ny,)c. (3.26) 

Here = and 

{Hxi ■ Uy-^Ux^ ■ 'l^y'ilc (3.27) 

1 1 

is the usual connected 4-point function, related to the previously used second if-moment 
as indicated. In terms of Wi and W2 the Binder cumulant reads 

^ ^+1 E.,,(w^±,i(^,i/) + i)]^ ■ ^' ' 

The Binder cumulant has accordingly a large N expansion of the form 

^± = E(]vTl)^^±-(^'^)^ (3.29) 

The coefficients can obviously be expanded in terms of the coefficients of the 2- and 
4-point functions summed over all arguments: 

^±,.(A,1^) := J][w^S(^,y) + A-^5,o] . (3.30) 



x,y 
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The two lowest orders 

U±fl = 2 - a^lw±fl , 

^±,1 = -o-±fo«^±,i + 2o-±,iCT^Jw±,o , (3.31) 
involve only functions already computed in the previous subsections. We have 

V 



w±,o = 2a^o 



2 



2V 



a;4n±(0) ' 



(3.32) 



and 



2V 



,1 = ~2g±^ — w±,o + 2 ^ r±{u,v)r±{w,z) 



(3.33) 



xA±(ii - w)[-D±(u - w)D^^(v - z)+ D±(u - v)A±{v - z)D±(w - z)] 
- ^ r±{u, vfA±{u - v) , 



where 



r±(x, y) := ^ W^},{z, w; x,y) = ^\l- U±{0r'D±,2{x - y) 



z,w 



(3.34) 
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4. TD limit of spin and current two-point functions 

Up to this point we have been considering both the compact and noncompact models 
in arbitrary dimensions d > 2. In the following we restrict attention to d = 2. Also 
numerous articles have dealt with the expansion of the compact model so we will 
in this and in the next section restrict attention to the noncompact case and drop the 
minus (— ) suffix on all functions. 

The results summarized in Section 2 seemingly suggest a simple relation between the 
compact and the noncompact models. It is important to stress, however, that the rela- 
tions hold only on a finite lattice and for invariant correlators. Physical quantities arising 
after taking the limit of infinite lattice size (thermodynamic limit) turn out to be very 
different in both systems. We will illustrate this fact by studying the thermodynamic 
(TD) limit of the coefficients in the expansions of the correlators computed in the 
last section in the noncompact model. The very existence of the limit is non-trivial in 
this case because a; — > as — > cxd, specifically 



This means a 'coordinated' limit of lattice sums of the form '^ki fc„ • • • ; ^n) 

has to be taken, where fy via u carries an explicit V^-dependence. The gap equation 
(13. 4p effectively acts as a subtle infrared regulator whose usefulness is underlined by the 
result summarized in Subsection 2.2. As mentioned, the sums associated with individual 
Feynman diagrams will typically diverge in the limit. The issue is whether the infrared 
divergences cancel in the W^^ and the quantities computed in terms of them. 

In this section we discuss the limit of the spin and current two-point functions; the limit 
of the Binder cumulant is computed in Section 5. 



4.1 TD limit of the spin two— point function 

In the leading order the 2-point function has an infinite volume limit 



where here and in the following means integration over the Brillouin zone J^^ (l^' 
The infinite volume lattice propagator V{x) is a remarkable function which has been 
discussed in detail by Shin [13]. At every lattice point it is given by an expression 




(4.1) 




(4.2) 
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of the form ri{x) + r2{x)/7i where are rational numbers. As \x\ 
logarithmically: 

V(x) ~ --^ 
47r 



(Inx^ + 27 + 31n2) + 0(|x| 



oo it diverges 



(4.3) 



where 7 ~ 0.577 is Euler's constant. 
The next-to-leading term is given by 



W^i'^(^'0) = :i5^(e^^^--l)W^«(p) 



(4.4) 



with 



{Ep + 00) 



(4.5) 



and g = g_ as in 03.91) . For p = the TD limit of Wi'^\p) does not exist (similarly to 
the situation for the leading order): W^^^'*(0) ~ yinVlnlnV", reflecting the fact that 
1^1(0:, 0) is an increasing function of distance 

For p ^ 0, however the limit exists. To see this we first note that 11 (p) has a TD limit 
for p 0. Indeed using one insertion of fl5.17b ) below and the gap equation one can 
rewrite Il{p) as 



U{p) 
J{p) 



(Ep + 2uj) 



J + Jip) , Py^O^ 



Ek + Ep-k — E, 



-y 

V ^ {Ek + uo){Ep.k + uj) 



(4.6) 



Throughout we often use the symbol J to denote lattice sums which give rise to conver- 
gent integrals over the Brillouin zone upon taking the infinite volume limit. The limit 
noo(p) of n(p) is then given by 



noo(p) 



1 



Epvij)) ' 



v{p) :-- 



A 



+ ^00 (p) 



(4.7) 



with 



Joo{p) 



Ek + -E'p-fc — Ep 
k EkEp-k 



(4.8) 



The properties of the function Joo(p) will be important later on; we mostly need: 



Joo(p)>0 forp^O, Joo((7r,7r)) = 



(4.9) 
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and the behavior for p ^ 



Joo{p) = -7^[lnp^ - 5ln2] + 0{pHnp^] 



(4.10) 



Eq. fl4.10p follows from [13], the positivity in (14 .Qp follows from Appendix A of 0. A 
direct way to see Joo(p) > is by performing one of the integrations explicitly. This 
leads to the integral representation 



2n 



-th- - th- + ^^(^ + ^) f 4sh^(^)-Pi 
^2 2 shashfo Ush2(^)+pi 



(4.11) 



where a, 6 > are determined by sh| — , ^ 
of the last term in the integrand one then has 



sin ^1, sh| 



sin ^^"2^ I • For the numerator 



4 sinh 



a + h 



. -, a . h . { a + h 

smh - smh - cosh 

2 2 V 2 



. fki-Pi\ (Pi 
sm,- sm(^— cos 



a h 
8 smh — smh - 
2 2 



cosh 



a 



e cos 



Pi 



(4.12) 



where e = ±1. For fixed the integrand of ( 14. lip therefore is a monotonically 

decreasing function of p2 for Q < p2 < n. By symmetry the same must hold with the 
roles of pi and p2 interchanged, so that Joo(p) > Jooiij^, tt)) = 0, for all p 7^ 0. 

Returning to (14. 5p we rewrite the sum as 

1 



-E 



l^^n(A;)(Ep_fc + cu) 



Mp) - 



An(p) 



(4.13) 



1 ^ 1 / 1 1 



U{p - k) Yi{p + A;) n(p) 



Then using the properties of the 11 function (14. 9 p and (14.1 OP . it follows that Jn has a 
finite TD limit which we denote by Jn,oo- In Section 5 we show that also the TD limit 
of q exists. 



goo = - / v{k) = Jn,oo(0) . 
'k 



Putting the results together one sees that the limit of (14.50 is 

2 



wi'LiP) 



Er. 



[3{p) + \ \{p)\ , p^Q. 



(4.14) 



(4.15) 



i{p) := — [Jn,oo(p) - ^n,oo(0)] . 
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From (14.131) one also gets the 'sum rule' 



I^E,wi%{p) = 2 + ^q^ , (4.16) 



while / W^i,!i(p) = = Wi^l^{x,x), as required. 



It is instructive to compare now the continuum (i.e. small ap behavior where a is the 
lattice spacing) of Wi^^{p) with its counterpart in the compact model. In Appendix B 
the small p asymptotics of j{p) is determined. In terms of the (non-universal) constants 
gsW the small p behavior comes out as 

EpWi^Lip) - - In (- HpVT)) + 2g,{X) + ^^f^^^j^'^ ^ + Q {Mp'mr) , (4-17) 
where 

T = 32exp l^y^ . (4.18) 

Eqs. (14.171) . (14.181) illustrate in particular the nonperturbative nature of the large N ex- 
pansion in the noncompact model, despite the fact that in infinite volume the expansion 
is effectively performed with respect to massless fields, as it is the case in perturbation 
theory. The infrared regulator uj{X, V) clearly works very differently from a constant 
'small mass' regulator. 

The subleading logarithmic terms in (14.171) are difficult to determine on the lattice. 
Using a continuum cutoff instead and cutoff-normalized continuum momenta 

^ E - / ^(^' - ' -] 3 - := ^ , (4.19) 

the continuum counterpart of (14.151) can be found in closed form: 

2T^/Pl\ 4:71 1 



pI \tJ X Ht/pI) 



\n{T/pl) Jo |1 - s| ln(T/p2) _ Ins ' 

where Li(x) = ds/\ns. The integral in the last term is singularity free as the diver- 
gence of the 1/1 1 — s\ factor at s = 1 is removed by the Ins, and InT/p^ > Ins holds 
over the entire range of the integration. The asymptotic expansion of (14.201) comes out 
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as 




c, 



-n 



(-l)"+ir(n + 2)[2 - ((-1)"+^ + l)an + 2)] . 



(4.21) 



We add some remarks. First, note that fl4.2ip contains factorially growing terms both 
with oscillating signs (Borel-summable) and with constant phase (non-Borel-summable). 
Second, (14.211) is an expansion in terms of the running coupling 



which can be re-expanded in terms of positive powers of the bare coupling A. Com- 
paring the result with its counterpart in the compact model one may verify that both 
perturbative expansions are related simply by flipping the sign of A; see [5] for a gen- 
eral proof of this "perturbative correspondence". Third, the expressions (14.171) . (I4.2ip 
do not suggest a nontrivial continuum limit reachable merely by a multiplicative field 
and a coupling renormalization. For example defining a renormalized coupling A^ by 
y\2g47r/A _ ^2g47r/Ar ^ ^Y\_e infinite cutoff limit can only be taken by allowing negative bare 
couplings. For all A > both A^ and a{p) vanish for A — > cxd. Finally, the expressions 
(14.201) . (I4.2ip can also directly be compared with their analogues in the compact model, 
see [H] for the former. In the compact model the bare mass gap itiq = A^e~^'^/^ enters 
and for ni^/p^ = 32e~^^^^/pl <^ 1 (i.e. small A at fixed pc) expressions for the self energy 
Ep^Wi^dPc) are obtained related to (I4.20p . (I4.2ip formally by flipping the sign of A. The 
attempt to take the sign flip beyond the asymptotic small A expansion, however, would 
on the bare level produce a "hyperviolet" mass scale A^e^'*'^''^, difficult to interpret. In 
contrast, in a lattice formulation the exact "large N correspondence" summarized in 
Subsection 2.2 exists. 



4.2 TD limit of the Noether current two— point function 

To obtain the infinite volume limit of the leading order contribution to the current 
correlation function (I3.20p we decompose it according to: 



(4.22) 



2n ■ HT/pl) 1 + AinA^ 



jj^Jiq) = 2exp f'-[q^ - q,]] [A^,{q) + sin(g^/2) sin(g,/2)n(g)] , (4.23) 
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with 



^ ^_ 1 -^ sin(p + g/2)^sin(p + g/2)^-sin(g^/2)sin(g^/2) 

Now we have seen in the last subsection that n(g), g 7^ 0, has a finite hmit, and so 
obviously does A^yi^q): 

^oo,i,u{(i) = I -7r4 — [sin(p + g/2)^ sin(p + g/2)^ - sin(g^/2) sin(g^/2)] , 

(4.25) 

which is independent of A. Thus the infinite volume limit of J^J{q) exists. 

In the next order we have the result ( 13.24p . We now show consecutively: (i) that X3 has 
a finite TD limit, and (ii) that Xi + X2 has a finite TD limit. 



(i) Noting the identity 



^ sin(p + g/2). 



p 

we can write X3 in fl3.25p as 

^ n = A V sin(pi + g/2)^sin(p2 + q/2)^Y{pi,p2) 

y2 (^^^ + ^) (^^^^^ + ^) (^^^ + ^) (^^^^^ + ^) ' I • J 

where 

y (pi, P2) = YTT^ T - TT7-T - W7-T ' (4-28) 

which vanishes when either = or ^2 = 0. Then we break up X3 

^3;fj.u = 77; ; [X^i^iu — X5;fiu " ^5■,u^l " ^Gi^iu] , (4.29) 

with 

^ .X _ 4 ^ sin(pi + g/2)^sin(p2 + g/2).F(pi,P2)(^g-^pJ(^c-^i ' 



P2y 



X. / N 4 sin(pi + g/2)^sin(p2 + q/2)^Y{pi,p2){Eg - EpJ 

Y (\ ^ 4 sin(pi + g/2)^sin(p2 + g/2)i.y(pi,P2) 
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Now the limit of X4 exists. Next for X5 we write 



^5;m.(5) = X7.,^M) + Xs;Aq) , (4.31) 

with 

4 ^ sin(pi + g/2)^sin(p2 + g/2). [^(^1,^2) - -q)] {E^ - E^,) 



P1,P2 



where (using the gap equation) 

r.i^l) - ^ E = J [1 + (4 + -)A-^] sin(g,/2) . (4.33) 

It is clear that both X7, and so also X^ have finite hmits. Finally we have 

Xe.^,{q) = Xg,^,{q) + X^o-Mq) + X^o-Mq) + 4i^(?> (lYMYM > (4-34) 

with 

4 X - 

X^;nv{q) = ^ 2^ sin(pi + g/2)^sin(p2 + g/2),. 

P1.P2 

^ [>"(pi,P2) - y(pi, -g) - y(P2, -g) + r(g, q)] 

(-E'pi+9 + CU) (-Ep2+9 + 



^ . ^ / ^ 4 V- sin(p + g/2)^[r (p, -q) - Y{q, q)] 

It follows that X^ has a limit, and the demonstration that has a finite TD limit is 
complete. 

(ii) We first rewrite the sum 

Xi;ti.v + X2-fj,v — —^q-A2-^u + Xii-f^^ , (4.36) 

with 

1 Ff,^{p,q) 



^ pT;^ n(pi - P2) (£^pi + w) {Ep,+uj) 
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Here 

p ^ sin(p + g/2)^sin(p + g/2)^ _ sin(g/2)^sin(g/2)^ 

To proceed we again write Xn as a sum of terms: 

Xii.,^.{q) = X,2-Mq) + Xis.^,{q) + SfiA^.^M , (4-39) 



where 



^l2■^l.Aq) - Tn. 2-^ — \ — ^ftt; — \ — : > 



and (/2 will appear in the next section): 

fs-=-Y . . / Ts . (4.41) 

We still need to consider the limits of Xi2;^,vi Xi^.^^, A2-^^u (actually the contribution in 
Xi + X2 involving A2 has a coefficient proportional to q — fi which vanishes in the TD 
limit). We note that these three functions can be written in the form T2-ni,[g]{q) where 
functions Ts-fj,^[g]{q) are defined by 

with g(j)) a regular periodic function with g{p) — g{—p) and finite at p = 0. Now for 
functions of the type 71 we have 

p ^ 

Er, + 2uj\V ^ E^ + uj 



1 9{p)Ff,u{p,q){Ep+q + u) 
Ep + uj 

p 



V 
p 
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which has a finite hmit since y Yl/p^tJ-viPiO) does, as can easily be seen using the gap 
equation. Next 



1 



{Ti-A9]{<i)+T:,,A9]{q)+T,,M{<i)} ■ 



Here 



Ts-Mgm -yi^ Ie;^^^ 



Ti-„A9]{<l) 



1 1 9{p)Ff,,{p. q) 



{Eg + uj)V ^ {Ep + co 



with 



Ffxuip, q) := {Eg + uj){Ep+q + a;)F^^(p, q) . 
Decomposing T3 further gives 

1 



Tz;A9M 



Ea+UJ 



{n-Mi^) + T6;,A9M} , 



with 



T - ^ 9{P)F^,U{P, q){Eg -Ep- Ep+g){Eg - Ep+g) 

T,;,A9\{q) -yl^ {Ep + uf 



1 1 ^ 9{p)F,u{p, q){Eg -Ep- Ep+g) 



E 



Eg+ivV^ {Ep + coY 



where T5 clearly has a TD limit. Now 

KAp, q) = F^.^p, q) + F+.^p, q) , 

with q) = ±F±{p, q): 

F-,i^u{P: q) ^ {Eg + uj) [sinp^cosp^cos(g/2)^sin(g/2)^ + ^ ^)\ 

-2 ^ sinpp cos qp sin(g/ 2)^ sin(g/ 2)y , 
p 

F+;u,u{P: q) = (-Bg + a;)sinp^sinp^cos(gf/2)^cos(gf/2)^ 



(4.44) 



(4.45) 
(4.46) 

(4.47) 



(4.48) 



(4.49) 



(4.50a) 



-^sm{q/2),sin{q/2)J{Eg + u;){pl+pl-lpy,) + 2Ep-Y,PX] ■ (^-^Ob) 
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So 



with 



t^-AqM = TT— - {T7-A9M + n-AgM + 5(o)r9;^.(g)} , (4.5i) 



^ P (^^ + ^) 

p p VP/ 

= -2 E cos g,- E ■ 

All of these quantities have a finite limit, in particular Tg.ni,{q) on account of the gap 
equation. A similar decomposition can be performed for T^-^uig), showing that it likewise 
has a TD hmit. So T2-^i,[g]{q) has a TD limit. It follows that -^i2;/ii/, -^i3;/ii/, ^2;/ii/ all 
have infinite volume limits. 

To summarize, wc have shown that the TD limit of the contribution to the two leading 
orders in the expansion of the current correlator jju) , s = 0, 1, exists. 



5. TD limit of the Binder cumulant 



In the following we evaluate the infinite volume limit of U in the 2-dimensional noncom- 
pact model to sub- leading order, i.e. the coefficients Uq, Ui in 

are computed for large volumes. Somewhat surprisingly we will find that the limit 
V ^ 00 exists and is independent of A! In Subsection 5.1 the large volume asymptotics 
will be evaluated analytically. As a test of the estimates and in order to have finite 
volume results to compare Monte-Carlo data with, we also directly evaluated the multiple 
lattice sums numerically up to L = 1024. The results are reported in Subsection 5.2. 
Finally, to preclude that the large N results are misleading we performed a Monte-Carlo 
study of [/(A, V) up to L = 384 for N 
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5.1 Analytical analysis of large V asymptotics 

Evaluation of the leading order coefficient is straightforward. From (13.311) and (13.3!^ 
one obtains 

Now the TD limit of ct;n(0) can be evaluated from 

- -n(o) 4 + i i: ^i-^ . -L V. + 1 + + o ( J^) . (5.3) 

where a(A) > 0. Note that (taking the A-derivative of the gap equation (13.41) ) this 
translates into a large volume asymptotics for a; = cc;_ of the form 



InV^ A In^V A2 V / J\n^V \\n^V 



From (15. 3p one has 

^7ra(A) 



87ra A / 1 
fo(A.V') = 2-l^ + 0(^ 



.In^l^. 

[/o(A,oo) = 2, (5.5) 

in stark contrast to the compact model (as discussed in the conclusions). 

The discussion of the sub-leading order is more involved and is best done in Fourier 
space. We prepare the following auxiliary functions 

1 v-^ 1 

V {Ek + ujY[Ep_k + 

Ils{p) = n,i(p), s>l, (5.6) 

where Il{p) := ^i{p) = A(p)~^ is the inverse of the Fourier transform of A(a;). Further 
we shall need in addition to (I4.4ip 



p,k 

Note 



/2i = n(0)g. (5.8) 
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The expression (13.291) for Ui we break up according to 

Uii = 2aiaQ^Wo, 



U- 



12 



-(To ^Wi 



[A + B + C + D]. 



Using di = 2Vu-^{-q + fi) and wo = 2Vuj~'^{V - LJ-^U{Oy^) one has 



11 



vu{o)J ■ 



In Ui2 the term A comes from the dw^/duj derivative in wi and is given by 



n2(o) 

VuoY[{<S) ' vn(o)^ 



+ 



The term B corresponds to the w!f'^ D /\D ^W^f* piece in wi and is given by 



B = -S[Vuj- ^ItttV/i - 4^/31 • 



c<;n(0) 



n(o)^ 



The term C arises from the w!f^ D ^^DW^' structure and reads 



^ n 



n2(A;) 



(A;)nn(o) ^fc + cjj 



Finally D comes from the Wf^mif^ term 



:D 



In the V ^ oo limit certain terms are superficially divergent, so it is best to 
terms where such divergences cancel. Specifically we rewrite A + B as 



1 



UJ 



2^2 (^ + ^) 



+ 



Vujm{o) n(o) 



(/in2(o)-/3i) 
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Further useful combinations to discuss the hmit turn out to be 

1 ^ 1 rn2(fc) 1 



q-fi 



(A;)Ln(0) (Ek + uj) 



/3i-n2(o)A=iEH 



n2(o) 

Eu + ui 



(5.16a) 
(5.16b) 
(5.16c) 



For the evaluation of the TD limit then mainly the combinations in square brackets 
in (15.161) have to be studied, for large volumes. A naive replacement of the lattice 
sums by integrals over the Brillouin zone with uj = would produce infrared divergent 
integrals. The strategy in the following will be to evaluate the large volume asymptotics 
of the functions lis, s = 1,2,3 and 1122 by repeated insertion of one of the following 
decompositions of unity 

1 = ^ [Ep - Ep_k + {Ep-k + uj)] , (5.17a) 
rjp + u} 



Ej, + 2uj 



[Ep -Ek- Ep^k + {Ep-k + uj) + {Ek + to)] 



(5.17b) 



until terms corresponding to infrared convergent integrals over the Brillouin zone arise. 
The volume and the momentum dependence of the additional pieces picked up in the 
process (which may diverge as — oo) can then be studied analytically. For Hi = 11 
itself only one insertion of fl5.17b ) was needed and no divergent piece arose, see ( 14.61) . 
Proceeding similarly we derive the relations 

= ^ [Xstip) + n(,_i)t(p) + n,(t_i)( 



{Ep + 2uj) 



sA>l 



Er, — Eu — E, 



p—k 



V ^ {Ek + ujyiEp-k + u;y ■ 

Applying it for the case s = 2,t = 1 (and noting Hsoip) = ns_i(0)) we get 

1 



Ep + 2uo 



[X2i(p) + n(p) + n 



This can be rewritten as 



1 



n(p) 



^2{P) 

Ln(o) 
1 



Ep + u 



n(o) 



X2l{p) 



Ep + 2uj {Ep + 2uj)Ii{p)\ n{p){Ep + uj){Ep + 2uj) 



(5.1^ 



(5.19) 



(5.20) 
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from which q — fi can be evaluated. In a first step one finds 



1 ^ 1 1 ^ Jsik) 



Here 



1 



y^n(A;)(Efc + cu)(i?fc + 2a;)' 

'^^^ ^ V^^U{k){E,_, + u;){Ek + 2u;)- ^^'^^^ 
The function Jslp) has a finite hmit which for small p behaves as 

J3,oo(p) = Ep[j3 + 0{1/W)] , 

h = / 4?(cos h - COS . (5.23) 

This can be used to show that 

n(0)cu/2 = -^Inlnl^ lnl^+ ^lnl^ + 0(lnlnl^), 

n(0)(g-/i) = -— lnlnFlny+^lny + 0(lnlny), (5.24) 
where the constants are related by 

?i + ?2 = ia ■ (5.25) 

Indeed, using the fact that E]^Y\.[k) has a finite limit which scales like — ^In/c^/T for 
A;^ — > 0, one readily verifies the first equation. Further 



n(o) "Vln^l^. 

+ vJuio) ^ ^^ + ^(l^)' 

y^(Efc + 2cu) ~ Fa; \\nV ) ' 

V^ji^. - J^.3.ny.0(l). (5.26, 



26 



Inserted into fl5.2ip gives 

n(0) [q-fi + uh] =j3^\nV + 0(ln In V) , (5.27) 
and hence the first equation in (15.241) . Combined with (15.31) we arrive at 

>-/')(l.l|.-0--"-o(w)' 

using 

Since Ui = Uu - {A + B + C + D)/{uj'^V'^) one concludes from flOSD and (15151) that 
the softly decaying Uu terms in Ui cancel. 

We proceed with the evaluation of C. To this end a more detailed evaluation of X21 
defined in (15.181) is needed. In a first step one obtains 



fi k 

+ { 2^ cos y + J2(P) + 2^COS YCOSP/,-2_^ 



1 (sin ki + sin 



2 



\2 ■ 



(Ek + iur 



k 



-64sin^^sin^^ W'^^^^]. (5.30) 



where 



J ( \ _ }_ ST^ (-Ep — Ek — Ep_k){E.p — Ep_kY 



-2< J2,oo(p)<0, J2,oo(p) = -2+j2p' + 0(/), j2~0.93. (5.31) 

As V — s> 00 the first curly bracket diverges logarithmically while the second one is 
convergent. Using 

1 ^ (sin^sinf)^ ^ J_^f]_ 
{Ek + ijf 327r \V 



1 (sin A;i + sinfc2)' 1,.. 1 1 

^ (Ek+uj)^ 47r ^ ^ 4 47r Vlnl^/ 



?2 
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with a(A) as in (15.31) one finds 



--7;^(,cospi - COSP2J 
p 



(5.33) 



+ 



^2,00 (p) + ^ "^{cospf, + cos2p^) + sin^P/. + O (j^) 



For small momenta this behaves like 

1 



X2iip) 



^ lnl^ + a(A) + - ^ 
Utt ^ ' 2tx\ L4 



1 /. 5 1 



P2 Pi 



+ o{pi,P2) 



(5.34) 



For the evaluation of C it is useful to rewrite (I5.18P for p 7^ as 



n2(p) 



n(p)Ln(o) Ep + cjJ (Ep + cu)n(o) 



where 



n(p) n(p) 



X2(p) := 1 



+ 2uj 



UJ 



X2{P) + 



_ Ujli2{p) X2l{p) 

n(p) n(p) 

X2i{p) Ep + uj 



U{p) Ep + 2uj 

(n(p) - n(o)) . 



(5.35a) 
(5.35b) 

(5.36) 



In the first term the fact that the numerator is proportional to Ep is important for the 
eventual decay properties of C. Further X2{p) is bounded by a function of order InV^ 
in the volume and it has at most logarithmic singularities for p — 0. The latter is 
consistent with 

ujU2{p) cun(O) 



U{p) EpU{p) 
where the 0(1/1^) piece comes from uX2i{p). 
Inserting (15.351) into (I5.13P gives 
8 



Oil/V), 



(5.37) 



C 



ujU2{0)\^ 8V 1 >^ 



m 



u{oyv^^{E, + ujr 



iUn2ik) , X2lik) 



+ 



.(5.38) 



Using the known behavior of the constituent functions for large V and small momenta, 
one can verify a decay of the fornix 



C 



O 



In^V^ 



O 



1 



(5.39) 



^Here and later on we indicate the form of the sub- leading term, without however (in a slight abuse 
of the O symbol) presupposing that its coefficient is nonzero. 
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We proceed with the D term, where 1122 (p) enters. A useful representation is 

+ (^^^^["'4(p) + 2J(p)-2A'i(p)], (5.40a) 

This can be used to determine the large V behavior of the term D in (15. 141) . We begin 
by separating the zero mode in the relevant combination 

1 f _ ( 1 ^ ^ ( ^^n22(o) 

n(o) - 2/2 + 1 - (^1 - y^2n(o)j " vcj2n(o) n(o^ 



1 1 



1 >^ 1 n22(A;) 2 

On account of 



-!(*:) L n(0) (Bi + uj) 



= 01 ' 



cj2-^n(0) Vln\/ 

cu'n22(o) ^ ^ c^n2(o) 
n(o) n(o) ' "Vin^v^ 



Vln^V^/ 



the zero mode pieces are 0(1/ In V^). For the last term on the right hand side of fl5.4ip 
we introduce the shorthand 81 + 82- Upon insertion of fl5.40p we write S'l for the part 
coming from the 11 (p) + 11(0) piece in (15.401) and 5*2 for the rest, 



^1 



1 ^ 1 2^1 



n(0) V^^{E, + 2uy V ^ n(A;) V {E, + 2ujf {E, + ujf 



0(^j. (5.42) 



The term 82 reads 



^ nk7g n(.)(E! + 2.)3 [^-(^) + ^^(^)-^^^(^)] 

+ 6 i>p X2m 

n(o)v^n(A:)(E, + 2c.)2' ^ • ^ 
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and is checked to behave as 



Together 



^ 0f-l-)+0f-i-). (5.45) 



It remains to consider A + B. In view of f l5.15p and fl5.28p we know 

SO that only /31 — n2(0)/i is still needed. Using (15.181) for the case s = 3, t = 1 we obtain 

/31 - n2(0)/i = ^— - ^^03(0) - uoli^i^)] + - V , 5^'^^^ r 

Vujm{Q)V ' '\ V j^^Tl{k){Ek + 2u) 

+ T^^n(fc)(E, + 2.;) ""^ ^^^^Vj^^Ii{k){E, + 2u){E, + uj)- ^^'^^^ 

Since 

a;2n3(0) - cjn2(0) = 0(1^2) ^ (5 48) 

the zero mode piece scales like 0(y/ Inl^). For the second term we observe 

with J3 as in (I5.22p . As a consequence this term in (15.471) scales like OiV) for large V . 
In the last two terms we insert (15.191) to get 



1 V 1 r. I ^2i(fc) i 1 1 

y^(E, + 2^)2[ ^ n(A;) J^F^n(A;)(E, + 2^)2 



n(0)-a;n2(0)^'= + 2^ 



Ek + uj 

(5.50) 

The very first term is OiV \n.V), the one involving X2i{p) is 0(1^/ In l^), and the last 
one is O(V^lnlnV^). Together 

/31 -n2(0)/i = 0(V^ln\/) + 0(rinlnV^). (5.51) 

For A + B this results in 

_i_„,,H,„,0(J^).0(!|^). (5.5, 
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Combining (15. 9p with (I5.52p . (15.391) and (15.451) we arrive at the conclusion: 

o(J^).o(!BiliK),o(^). 

f/i(A,oo) = 0. (5.53) 
For the TD hmit of the full Binder cumulant the result (I5.53P amounts to 

U(X,oo) = ^— + o( ^ ^ ^ . (5.54) 
^ ' ' N+1 V(A^+1)^/ ^ ' 

This result will be backed by Monte-Carlo simulations in Subsection 5.2. Potential impli- 
cations for "criticality" and "triviality" of the theories are discussed in the conclusions. 



5.2 Direct evaluation of lattice sums 



Both as a check on the previous analysis and in order to have finite volume data to 
compare Monte-Carlo data with, we also evaluated the lattice sums defining Ui and 
several other quantities numerically up to L = 1024. Since O(L^) terms have to be 
summed and both very small (e.g. uj) and very large numbers (e.g. n(0)) enter high 
precision is needed. The summations were performed to 96 bit (26 significant figures) 
accuracy using the publicly available arbitrary precision MPFR library (www.mpfr.org) 
and for moderate L also with Mathematica. 

The results were found to vary with A such that for smaller A the presumed large V 
asymptotics sets in later. Below we present the results for A = 3; qualitatively those for 
other A values are similar. Due to the predicted occurrence of very slowly varying terms 
(e.g. of In In V^/ In type) one cannot expect that the genuine large V asymptotics can 
be unambiguously probed by direct summation. Nevertheless two or three parameter 
fits of the L < 1024 sums to the expected decay form are generally convincing. Table 
1 summarizes results for u, n(0) and some slowly varying quantities entering Uu. Here 
Uii is defined in Eq. (15.101) . q — fi is evaluated directly from (13.91) . (I4.4ip and from 
(15.16b ). /2 is defined in (I4.4ip . For example the leading asymptotics ull{0) ~ —ji^liaV 
and the coefficients in 

-/2~ilnln\/~-(g-/i), (5.55) 
z ui 

come out well in fits to the data. 



Table 2 presents results for the terms used in the breakup of f/i, see Eq. (15.91) . and the 
final result for Ui. The column for A + B again illustrates the need for high precision, 
as individually A and B are 2 — 5 orders of magnitudes larger that their sum. It also 
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L 




iq- fi ) luj 


J I 


10-^11(0) 


64 


-233.007495171 


-0.25538202207 


0.42884117188 


0.00451330538 


128 


-52.7199137013 


-0.23716419150 


0.37935605329 


0.02202558370 


256 


-12.0362079518 


-0.21754473457 


0.33427610710 


0.10558882967 


512 


-2.76867400687 


-0.19723910189 


0.29289048252 


0.49868470422 


768 


-1.17557207628 


-0.18522986021 


0.27016715051 


1.22915930848 


1024 


-0.64094863000 


-0.17669596989 


0.25464412160 


2.32558377444 



Table 1: Quantities entering \J\\ for A = 3; all given digits are significant. 



highlights that the analytical evaluation of the large V asymptotics is crucial. Even at 
L = 1024 the normalized A + B contribution is still increasing. On account of f l5.52p 
the decay of the combination {A + B)/{uV)^ — Uu should be faster. Indeed these data 
have a maximum at around L = 300 and then decay monotonically in a way fitted well 
by the predicted functional form. 



L 


lO^t/n 


103 {A + B)/{Vujf 


103 C/{Vujf 


10'^D/{Vujf 


Ui 


64 


-0.748704 


-4.65746 


9.13703 


0.387281 


-0.012353890 


128 


-0.565668 


-2.24328 


7.25786 


0.275149 


-0.010946420 


256 


-0.431162 


-0.73430 


5.92122 


0.189879 


-0.009688426 


512 


-0.330143 


0.23478 


4.92648 


0.132507 


-0.008595209 


768 


-0.282636 


0.63551 


4.45707 


0.108460 


-0.008027416 


1024 


-0.253114 


0.86471 


4.16358 


0.094567 


-0.007654026 



Table 2: Quantities contributing to C/i for A = 3; all given digits are significant. 



Finally we present a fit of the Ui data to the predicted decay form in fl5.53p . 



5.3 MC results for U 

Since the large N expansion is only an asymptotic expansion the higher order coeffi- 
cients in fl3.29p are not bound to be small, even in finite volume. At any given the 
truncated series could in principle misrepresent the exact U{X,V). In order to preclude 
this possibility we estimated U{X, V) via Monte-Carlo simulations. 

We have chosen to simulate a S0(1, 8) theory at A = 3 on lattices of linear dimensions 
L = 32, 64, 128, 256, 384. The simulations were performed in a fixed spin gauge (the spin 
at the origin was held fixed). The variable spins were updated by a Metropolis procedure 
tuned to achieve a roughly 50% acceptance rate. Equilibration and autocorrelation 
times for various observables have an enormous range: non-gauge-invariant observables 
in particular (e.g. (n^)) require extremely long runs, and on larger lattices fail to reach 
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Figure 1: Ui vs V. Fit to 0{\n\nV/\n^V) + 0{l/\nV) + 0{l/\ii^V). 

equilibrium even after billions of Monte-Carlo sweeps. The situation is much better for 
gauge- invariant observables, such as S-E entering Ui, see (13.261) . The fluctuations in 
this latter quantity determine the Binder cumulant, and are typically stable after a few 
million sweeps. Results for the quantity — U versus lattice size are shown in Fig. 2. 
This quantity is not monotonic, but reaches a maximum near L = 64 and then decreases 
quite rapidly. The decrease appears to be faster than the log-type decay found for Ui in 
Subsection 5.1, suggesting that the termwise large V asymptotics of the large N series 
(when formally treated as convergent) sums to a power-like large volume decay. 




Figure 2: Monte-Carlo results for 



N+l 



[/ VS y 



In summary, the numerical evidence suggests that the large N contributions to the 
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Binder cumulant beyond leading order (i.e. f/o(A, V)) may indeed vanish in the thermo- 
dynamic hmit. 



6. Conclusions 

Noncompact S0(1,A^) sigma-models are expected to be massless in contrast to their 
compact counterparts. The infrared problem therefore is nontrivial, especially in dimen- 
sion d = 2, and the goal of the present paper has been to gain computational control 
over the limit of vanishing infrared regulator. The large N expansion is well suited for 
this; in a lattice formulation the dynamically generated gap is negative and serves as 
a coupling dependent infrared regulator which vanishes in the limit of infinite lattice 
size. The cancellation of infrared divergences has been demonstrated in c? = 2 by ex- 
plicit computation of a number of physically interesting quantities defined in terms of 
invariant correlation functions: the spin and current two-point functions as well as the 
Binder cumulant, all to next to leading order. A complementary result is [15] where 
a noninvariant observable was shown to have a finite thermodynamic limit in d > 3 
beyond large N. In d = 2 we expect that a 'large A^' counterpart of David's theorem 
[7] can be established, showing that infrared divergences cancel termwise in the large N 
expansion of invariant correlation functions to all orders. 

To discuss our result for the Binder cumulant let us first recall the situation in the 
compact model. In the notation of Subsection 3.4 one has there 1/11+ (0) ~ inu^, 
in the thermodynamic (and continuum) limit, so that VU = 87r/[(A^+ l)ct;+]. Taking 
= 1/ ^/u^^ as the definition of the correlation length, this gives the familiar result for 
the renormalized coupling g^ = VU/^"^ = 8n/{N+l), to leading order. See also ref. [3] 
for a direct continuum computation to sub-leading order, with the result (A^-l-l)gr = 
87r[l - 0.602033/(A^+l) + 0{1/{N+1)% 

In the noncompact model the zero momentum limits of invariant correlation functions 
are expected to diverge in the thermodynamic limit. Indeed, mostly this refiects the 
fact that they are increasing functions of the lattice distance (recall Ux-Uy > 1, always). 
Our results of Section 5 suggest however that the ratio entering U is finite, independent 
of A, and very close to 2/(A^-|-l). 

One can view this result as a manifestation of a "concentration of measure" phenomenon. 
For the ID lattice model with L sites it was shown in [16j that the functional measure 
has support mostly on configurations boosted by an amount increasing at least powerlike 
with L. In the thermodynamic limit the measure (or mean) is therefore concentrated 'at 
infinity', i.e. in the disc model of the hyperbolic geometry at the boundary of the disc. 
Though not proven in dimensions li > 1 it is very plausible that a similar concentration 
phenomenon will hold for the c?-dimensional functional measures. Indeed our result 
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on the Binder cumulant can be put into this context: First note that in terms of the 
normahzed average spins a" := S"/ ■ S, with = Xlx^"' o'^^ 

Var(a2) := {{a ■ a - {a ■ a)f) = -U>0. (6.1) 

In the thermodynamic hmit Var((j^) has been argued to vanish, which is natural if the 
components of a" are typically very large rendering the relative fluctuations ensuring 
(cr-a) = 1 and ((a-cr)^) ~ 1, negligible. Indeed to leading order of the large N expansion 
one finds (cr°) ~ y/ln V. Note that the constant 2/ (A^+1) can be interpreted as the value 
of t/ in a constant configuration and that the indefinite dot product is crucial here. 

Alternatively 2/{N + 1) — f/ is given by the ratio of the susceptibilities defined from the 
partially connected 4-point and 2-point functions W2 and Wi, respectively. Both diverge 
as V ^ 00, but the ratios entering the large expansion of U (viz wi/a^, 102/(7^, see 
Eqs. (13.281) - (I3.3ip ) vanish in the thermodynamic limit. This is compatible with a 
genuine factorization of the 4-point function but does not entail it. (Since the connected 
four point function {rix^ ■ ny^n^^ ■ ny^) ^ entering (13.261) does not take into account the 
nonzero one-point functions, the fact that it must be non-zero is only indirectly relevant 
for this.) 

Concerning local quantities, the analysis of the TD limit for the subleading term of the 
spin- two point function in Subsection 4.1 does not suggest the existence of a nontrivial 
limit as the UV cutoff is removed. Positive bare couplings are required for the large 
series to be an asymptotic expansion, in which case a naturally defined renormalized 
coupling vanishes as the UV cutoff is removed. The situation should be similar for the 
two-point function of the Noether current. 

Together, our results may be taken as an indication for "triviality" of the theory in the 
sector comprising S0(1,A^) invariant observables. If corroborated beyond the large A^ 
expansion this would be of significance in a number of other contexts, e.g. for a class of 
Kaluza-Klein theories or for the widely studied systems with AdSs x target spaces. 
The focus on invariant observables is certainly natural from the viewpoint of the compact 
models. In the context of the Osterwalder-Schrader reconstruction H], however, 
invariant correlators are not ideal for the non-compact systems, and the situation may 
well be be different when noninvariant observables are considered. 

Acknowledgements: We wish to thank E. Seller for many discussions, and M. Liischer 
for correspondence. The research of A.D. is supported in part by the National Science 
Foundation under grant PIIY-0554660. 
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Appendix A: Leading and next-to-leading order SD equations 



Here we tabulate the first few of the hierarchy of Schwinger-Dyson equations for the 
large coefficients r + s > 1. They can be obtained e.g. by first converting (13. ip 

into a system of equations for the exact Wr and then inserting the large ansatz (12.61) . 
The ensued recursive structure is summarized in Fig. 1. 

The leading 2-point function ( 13. 2p satisfies 



(0), 



\z=x 1 ^xy ■ 



(A.l) 



In the next order we have 

± [A^-Ui] Wi%,y) - XD±{x,y)A,wi%,z)U=, 
XWi%,y) + l\ ±XA,wi%,z;z,y)l 



6 



(A.2) 



where we have used the solution (13. 2p to the leading order equation (lA.ll) to simplify 
some terms. 



We see that to solve (\A.2\} we first need to solve the equation for the leading order 
4-point function: 



I z=xi 



±{A^, - uj±)W^l{xi,yi;x2,y2) - XW^l{xi - yi)A^^W^l{xi, z; X2,y2) 
= =F [6^^^^D±{yi - y2) + 5^^y^D±{yi - X2)] 

+A [5x1X2 + ^xry2] D±{xi - yi)D±{x2 - 2/2) ■ (A.3) 



The solution for is given in (13.51) . from which can verify that W)^'i in (13.81) solves 
(lAj). 

In the next order the equation for the 4-point function is 

± [Axi - uj±] W''^\{xi,yi;x2,y2) - XD±{xi - yi)Aa,^W^ll{xi,z;x2,y2)\z=x^ 



(1) 



= ±XA^^W^l{xi, u; V, yi, X2, y2)\u=v=xi 
±Aiyi°J(xi, yi;x2, y2)A^^W^l{xi, z%=^^ 

±AI^£i(Xi, 2/i)Ax,W^i°) (Xi, Z- X2, y2)\z=xr 

+A (5x1X2 + ^xij/a) |vr±°2(a;i, 2/1; a;2, 2/2) 

±/^±(x2 -2/2)H^S(xi,yi) ±Z^±(xi -2/i)iy£l(x2,2/2)} . 



(A.4) 
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One sees the pattern summarized in Fig. 1 emerging, in that the solution of flA.4l) requires 
knowledge of the leading order 6-point function. The latter satisfies the equation: 

± (A^.^ - uj±) W^lixi, yi, X2, 1/2; 3^3, 2/3) 

= ±W^l{xi,yi; X2, y2)A^,W^l{xi, z; X3, y3)\z=xi 
± W^l{xi,yi; X3, ys)A^^W^l{xi, z; Xa, y2)\z=x^ 

-S^r^^^^W^liyi, y2] X3, 2/3) - 5x1X3 (1/1 5 1/3; X2, 1/2) 

-^x^y2W^l{yi, X2; X3, ys) - S:,,y,W^l{yi, X3; X2, 1/2) 



(A.5) 



±A [4iX2 + 5^12^2] D±{xi - yi)w!fl{x2, 2/2; X3, ys) + D±{x2 - y2)W^'2{xi,yi] X3, ys) 



(0), 



±A [4ix3 + ^is/s] D±{xi - yi)W^Ux2, 2/2; X3, ys) + £'±(^3 - ?/3)lV]!'^(xi, yi, X2, 2/2) 



(0), 



where all the functions on the rhs are known from solutions of f lA.l|) and (lA.Sp . The 
solution is simply given by 

W^l{xi,yi;x2,y2;x3,y3) = ± ^ ^ 

tOi,W2,-!iJ3 21,22,23 



xWlJ^ixi^yi, wi, j(x2, 2/2; W2, ^2) W'±,2(a;3, Vs] W3, Z3) 

xD^\wi - ^2)^±^(U^2 - Z3)D^\W3 - Zi) . 

Using flA.6P one can verify that as given in (13.111) solves flA.4p . 



(A.6) 
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Appendix B: Continuum limit behavior of j{p) 



In this appendix we consider the small p behavior of the function j{p) entering the spin 
two-point function f l4.15p to subleading order. We have 

Sm J2;^ 



with 



jiip) = - [ -J—[v{k)-v{p)] , 

Jk J^k-p 

J2Ap) = [ ^ Hp - fc) - viP + k)] , (B.2) 
Jk 



h-Ap) = -7 Hp - k) + v{p + k)] 

4 Jk ^k 



where v(p) is defined in (14.71) and A;^ = 2 sin as usual. Note first j^ Ap) non-singular 
at J9 = 0: 

j3;,.(0) = -^Jnoo(0). (B.3) 



Next 



sin(A; -p)^,. 



J2Ap) = -2/ ^ -^'^ [y(k)~vip)] 

Jk J^k-p 

= - sinp^[2ji(p) + jiAp)] - 2 cosp^ j5;/.(p) , (B.4) 



with 



f k"^ 

jiAp) = -^Hk)-vip)] 

Jk ^k-p 



J5;m(p) = -ET^l'^ik) -v{p)] = 2y2smp^je.f,^{p) , (B.5) 



j.Mp) = m^^Hk)-v{p)] 

Jk ^k-p^k+p 



Noting 



j4;^(0) = -^Jnoo(0), (B.6) 



we have for small p^: 



J{p) - -3i{p) - 4^^^^^|?^j6;,.b) + ^ Jnoo(O) + 0{llW) . (B.7) 

hp 4 
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From (14.101) we have 



v{p) = a{p) + 0{p^) , a{p) :-- 



2tt 



\n{pyT) 



(B.8) 



with T defined in f l4.18p . We now consider two corresponding integrals 

- 9{c^ - e) ^^^^^ _ ^^^^^ ^ 



jiip) 



{k -pf 



k {k-py{k + pf 



[a{k) - a{p)] 



(B.9) 



where denotes d^A;/(27r)^ and c is a momentum cutoff T > > p^ . These give 
the leading small p^ contribution because 



3i{p)-3i{p) = t^i(c) + 0(l/lnp2) , 



with 



iP) ip) 



fi(c) 



V2 



1a 

'2 



Vl{c) + -V2 



0(1/ In/) 



(B.IO) 



A;2 



First using 



. El 



2n 



2n 



(t + COS0) 



> 1 



(B.ll) 



(B.12) 



we can do the angular integrations in the j functions, setting without loss of generality 
p ^ (p, 0): 



ji{p) -- 

J6;00(P) = ■ 

Noting > p'^ we obtain 



da; 



2 Jq |x — p^I 
dx 



V Jo 



ln(x/T) In(pV^) 



X + p 



\x — p^ 



ln(x/T) ln(p7r) 



(B.13) 



ji{p) 



21n(p2/T) 



[S,{p\T) + S2{p\T)+S,{p\c\T 



(B.14) 
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with 

Iny 



S,{p',T) = f dy 
Jo 



^ (1 - y) In(pVr) ' 



'1 



|/ln[(l + |/)pVT] 



Now for small p 



2 



oo 



'"(^Vr)^ ' [ln(j>^/T)l» • 

^^(^^^) - i^E4"'p„. (B.ia) 



with 



Jo y 

• • (0) 7r2 (0) 7r2 , 

giving s\' ^ -\ ,sy ^ Next 

53(p^c^T) = 54(p^c^^)+55(p^c^T), c2<r. (b.is) 

Here 

' ' ' ^ " ii ^1 + ln[(l + y)pVT] ~ L z Hzp^/T) 



and 



ln(l + i/) 



1/(1 + 1/) ln[(l + y)pVT] 
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with 

-(-1) ^y^iYT^^ (^-'^^ 

giving sP = fj + |(ln2)2 .... So by flBlOl) and flRlil) 

Jib) ~ In (- ln(2pVT)) +^^1 + ([ln(pVT)]-i) 

^1 = t;i(c) + ^ln(-ln(cVT))+i, (B.22) 
which is independent of c. Similarly 

Mooip) = J^^^ [^6(p', T) + ^2(p', T) + Ss{p\ c\ T)] , (B.23) 

with 

= i (1 - !/) l'n(p%/r) ~ . (B.24) 



where 



(0) 1 

giving ' = 1 - -g- , . . . . 

Putting all the results together we obtain (14.171) with 



g2 = 9i + \v2 + \jiioo{^). (B.26) 
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Appendix C: Large with two auxiliary fields 



The results for the large expanded correlation functions in the noncompact model have 
in Section 3 been obtained via the large correspondence summarized in Subsection 2.2. 
Direct large computations in the noncompact model can be based on the following 
generating functional ^ 

e^pW^lH] = exp|-^^//,,}Ar/"j]d«,exp|-(Ar+l)5_[«,i/] 

x,y x^XQ 

S-[<^,H] = lTrlnA + ^^a.-i-(I-X'.o' 

x^xo 



■^xy ^xy ~l~ '^i^Olx^xy ~l~ ^ Hxy -\- HX5xy5xxQC>lxQ • (^"1) 



Here A is the matrix obtained by deleting the xo-th row and column of A or A. The 
formal expansion based on (IC.ip is not a valid saddle point expansion but it does produce 
the correct expansion coefficients and is related to its counterpart Wl [H] in the compact 
model by the involution ax i— > — a^, A i-^ —A. The functional (IC.ll) thus provides a simple 
heuristic way to understand the large correspondence. In contrast to the compact 
model, however, is not equivalent to the original generating functional 

Here we outline how fIC.ll) can formally be obtained from the formulation of the large A^ 
expansion with two auxiliary fields introduced in [4j. We begin by dualizing the 'spatial' 
spin components Ux-, x G A, as one would do in the compact model. Indeed, the A^ 
spatial components n^;, x ^ xq, enter (12. 4p with the 'good' sign; their 'dualization' gives 

expiy_[i7]= exp {-Ij^ j \{dnl5{nl-l) (C.2) 

x,y X 

X / JJda,exp|-^Trlnl-i(Ar+l)^a,|exp| + ^^^^^n°I,j,nj}. 

x^XQ x^xo x,y 



Here A arises due to the constrained Gaussian integration. Note the small but crucial 
differences to the compact model [5]: only A^ copies of TrlnA occur so far and the sign 
of the '^x^xo term is fiipped, as is the sign of the Hxy term in Axy. Most importantly 
the kinetic term in the last exponential has the wrong sign, which is why one cannot 
naively interchange the order of the integrations over ?t,°. and ax, x xq. To proceed we 
assume that in a large A^ expansion the replacement 

/ ll^<6{nl-l) = {-r\/' J l[dVxS{Vxo), n: = nx + zvx, (C.3) 
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is legitimate, for certain "saddle point" configurations fix, with = 1. With this 
replacement the kinetic term acquires the good sign. After the additional re-routing 
= —ioJxl (2A) + the saddle point conditions dS/drjx = dS/d^x = lead to 

nl = l + XDxx , (An)^ = UxUx , x 7^ Xq • (C.4) 

(See Eq. (5.18) of [1], with Ux = '2\dx, and correcting the sign in the second formula). 
Here Dxy = {M~^)xy with Mxy := —Axy + SxyUJx- Since —A is a positive operator 
it follows from the second equation in (]C.4p that the position dependent Ux must be 



UxUl. 



predominantly negative: < — ^^^^.(An)^ = — 
To leading order the spin two-point function is given by [1] 

{Hx ■ Uy) f,s. = rixfly - \Dx,y , 



D.,y ■■= Dx,y - , (C.5) 



XO,XO 



where we write momentarily ( ) j g for the average computed with the fixed spin measure 
in fl2.4p . The quantity fix then is the nonzero expectation value (^°)/.s. to leading order 
in 1/(A^ + 1). 

To make contact to the gap equation (13. 4p in Subsection 3.1 we now first replace flC.4p 
by a simpler gap equation with constant u and n, 

n'-\D'{0) = l, n'u; = --, with D' (x) := - ^^—^ ' (^.6) 

These are the saddle point conditions arising from a translation invariant gauge fixing 
of the functional integral (see Eq. (5.6) of [1]) and imply — AD(0) = 1, in accordance 
with (13. 4p . Given a solution uj^fi of (1C.6P we claim that 

fix := -XD{x - Xo) , iOx := uj + X5x,xq , (C.7) 

is a solution of (]C.4p . The equation {Afi)x = {uj + \6x,xo)nx is checked using —\D{0) = 1. 
To verify the first equation in (1C.6P it suffices to observe that 

Dx,y = D{x -y)- D{x - Xo)D{y - Xo)/D{0) . (C.8) 

This canjDe seen as follows: suppose that invertible matrices M and M axe related by 
Mxy = Mxy — cSxySxox- Theu the inverse of M is related to the inverse of M by 

{M-%y = {M-%y + iM-%xoiM~')yx, - (C.9) 

1 - c{M-^)xoxo 
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Applied to M^y = -A^y + UxSx y = M^y + XS^yS^^xo, for Ux = uj + X5x,xo, without yet 
assuming the gap equation flC.6p . this gives first 



Dxy = D{x -y)- Y^Xd(0) - ^o)^(^ - ^o) , (C.IO) 



and then fIC.SP from fIC.Sp . In particular no pole occurs for the quantities in fIC.SP at 



AL)(0) = -1. Using dOSD and AL'(O) = -1, it follows 1 + XDx^x = X^D{x - x^f = nl, 
while the sign in ( ]C.7I) is fixed by —XD{x) > 1. 

Inserting (1C.7P into (IC.Sp we arrive at 



{ux ■ ny)i,s. = -XD{x - y) = n - XD\x - y) = {ux ■ Uy) tr^^ns , (Cll) 

where the right hand side coincides with the spin two-point function computed to leading 
order in 1/(A^-|- 1) in the translation invariant gauge |1] and with that of Subsection 3.1. 

In summary, the leading order results with two auxiliary fields and the two gauge fixings 
considered (fixed spin and translation invariant gauge) are related via ( 1C.7I) . Both saddle 
point equations ( 1C.4I) and ( 1C.6I) imply the version — AD(0) = 1 used here, but in addition 
provide the interpretation of fix = {n'^)f.s. and n = (n°)trans, as the averages of n° with 
respect to the respective gauge-fixed functional measures. Note that fix approaches 
as |x — Xq I becomes large and that Yl^Ux = Vv?. Although the invariant two-point 
functions coincide to leading order in the two gauges, the results for the noninvariant 
quantity are very different, (n^^f.s = {nl)l^^. 

Equipped with this interpretation of fix we return to ( 1C.2I) . Subject to the assumption 
(1C.3I1 one can proceed by interchanging the order of integrations, which results in the 
Gaussian 



^xy^y 



X x,y X y 

= Const (det A)-'/' exp { [< (A-')-^^^ - n,.I,.,n,] } . (C.12) 

x,y 

Using also Uxq = 1 and substituting back into (1C.2P we arrive at ( IC.ip . 
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